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Abstract 

The electronic structure of SrBi2Ta20g is investigated from first-principles, 
within the local density approximation, using the full-potential linearized aug- 
mented plane wave (LAPW) method. The results show that, besides the large 
Ta(5d)-0(2p) hybridization which is a common feature of the ferroelectric 
perovskites, there is an important hybridization between bismuth and oxygen 
states. The underlying static potential for the ferroelectric distortion and the 
primary source for ferroelectricity is investigated by a lattice-dynamics study 
using the frozen-phonon approach. 

Pacs numbers: 77.80.-e, 77.84.Dy 
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I. INTRODUCTION 



Ferroelectric materials can display a wide range of dielectric, ferroelectric, piezoelectric, 
electrostrictive and pyroelectric properties. The potential utilization of these properties in a 
new generation of devices has motivated intensive studies. For example, the high dielectric 
permittivities of perovskite-type materials, like Sr x Bai_ x Ti0 3 , can be advantageously used 
in dynamic random access memories (DRAM), while the large values of switchable remanent 
polarization of ferroelectric materials are suitable for non-volatile ferroelectric random access 
memories (NVFRAM) ||-||. 

The most popular ferroelectric materials for nonvolatile memory applications are 
PhZr^Tii^Oa (PZT), because they have a high Curie temperature and large remanent 
polarization. However, these materials have serious fatigue degradation problems which can 
be solved by the modification of the electrode. An alternative approach to control the fa- 
tigue problem in ferroelectric capacitors is to use other ferroelectric materials. In recent 
years, SrBi2Ta20g (SBT) has emerged as an important candidate for non-volatile ferroelec- 
tric memories |5||J. It exhibits many desirable properties in order to be considered as an 
important component of the memory devices under development: almost no fatigue after 
10 12 switching cycles, good retention characteristics, low switching fields and low leakage 
currents. 

Bismuth-containing layered perovskites have been found to be ferroelectric by Smolenskii, 
Isupov and Agranovskaya [[jj . These materials belong to the family of Aurivillius compounds 
with a general formula (Bi 2 02) 2+ (A m _iB m 03 m+ i) 2 ~, consisting of m perovskite units sand- 
wiched between bismuth oxide layers || (here A and B are the two types of cations that 
enter the perovskite unit). 

It is well known that the Aurivillius composition SrBi2Ta20g (SBT) has ferroelectric be- 
havior at room temperature. Its crystal structure was investigated by Rae and co-workers || . 
The room-temperature structure is orthorhombic (space group A2xam) and the primitive 
cell contains 28 atoms. The lattice parameters of the conventional unit cell are a = 5.531A, 
b = 5.534A and c = 24.984A. As a result, its electrical properties are expected to exhibit a 
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high degree of anisotropy. For instance, it was observed that the ferroelectricity along the 
c axis is absent or very low [Hj, with a being the polar axis. More recently, this property 



was confirmed at the submiscropic level by piezoresponse scanning force microscopy [11 . 

The complex crystal structure of SBT can be described in terms of relatively small per- 
turbations from a high-symmetry body-centered tetragonal structure (space group symmetry 
I4/mmm, a = b ~ = 3.9lA), which contains only one formula unit per primitive cell (as 
in other orthorhombic structures, the a and b axes are rotated by 45° with respect to the 
tetragonal case). 

Although the orthorhombic symmetry is responsible for the ferroelectricity in SBT, the 
parent tetragonal structure shown in Figure 1 provides a convenient simplification for visual- 
izing and dealing with the complex SBT structure. As a matter of fact, most bismuth layer 
compounds have been reported to be pseudotetragonal, with tetragonal symmetry above the 
Curie point and orthorhombic symmetry below it (for SBT the Curie temperature is 608K 
). So the parent tetragonal structure of SBT could be related to the crystal structure of its 
paraelectric phase. Two main distortions from the tetragonal prototype structure lead to the 
orthorhombic structure. First, the ions displace along the orthorhombic a axis ([110] axis 
of the tetragonal structure). Second, the Ta06 octahedra rotate around the a and c axes. 
The first factor is directly responsible for the observed macroscopic spontaneous polarization 
along the a direction. 

In spite of its technological importance, theoretical studies on SBT are limited due to the 
complexity of its crystal structure. The band structure of SBT was initially calculated by 



Robertson et al. using the tight-binding method [12], where a highly simplified orbital basis 
set was used to reproduce the main features of the bonding. They used an orthogonal basis 
of O p, Ta d, and Bi s and p orbitals, and no orbitals on the Sr. The interaction parameters 
were found by transferring them from established band structures of other compounds. 
Their results show that both the valence and conduction band extrema are composed of 
states localized mainly on the Bi-0 layer. It was argued that the Bi202 layer dominates the 
electronic response (band gap, effective masses, etc.), while the ferroelectric response largely 
originates from the SrTa207 perovskite blocks. 



Several other speculations have been made about the origin of ferro electricity in SBT 
and related compounds. Among them are the central role of Sr +2 ion diplacements [13], the 
off-center position of Ta +5 ion relative to its octahedron of surrounding oxygens \T%\, and 
the movement of the Ta-0 plane relative to the Bi-0 plane ||. 

More recently, the difference between the electronic structures SrBi2Ta20g and 
SrBi 2 Nb20g (SBN) was investigated using the discrete variational X-a cluster method |Tp| . 
The results indicate that the difference in the remanent polarizations in SBT and SBN is 
due to the different displacements of Ta and Nb at the B sites of a pseudo-perovskite layer, 
and not due to differences in the displacement of the other ions. However, one cannot ex- 
pect high-accuracy results for relaxations and electronic states in crystalline materials from 
cluster methods. Thus there is a need for the application of highly precise self-consistent 
band structure methods. 

In this work we present a first-principles study of the electronic structure and ferroelec- 
tric instability in SBT. The calculations were performed within the local density approxima- 
tion to density functional theory, using the full-potential linearized augmented plane-wave 
method. 



II. METHOD 



The calculations presented in this work were performed using the full-potential linearized 



augmented plane- wave method (LAPW) method (see, e.g. Ref. 0]) with the addition of 
local-orbital basis functions [|17] as implemented in the WIEN97 code |[L8| . Exchange and 
correlation effects were treated within the local density approximation (LDA), using the 



parametrization by Perdew and Wang [19 



The muffin-tin sphere radii Ri = 2.0, 1.8, 2.3 and 1.5 a.u. were used for Sr, Ta, Bi and 
O, respectively. The value of the parameter RK max , which controls the size of the basis set 
for the wavefunctions, was chosen to be 7.3 for all the calculations. This resulted in well 
converged basis sets consisting of approximately 2500 LAPW functions. For the Sr-4s and 
4p, Ta-5s, 5p and 4/, Bi-6s and 5d, and 0-2s states local orbitals were chosen in addition. 
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Integrations in reciprocal space were performed using the tetrahedron method. We used a 
6x6x6 mesh which represents 28 k-points in the irreducible wedge for the body-centered 
tetragonal structure. Convergence tests indicate that only small changes result from going 
to a denser k-mesh or to a larger value of RK max . 



III. RESULTS AND DISCUSSION 



A. Electronic structure 



The electronic structure of SBT is calculated for the tetragonal parent structure. We 
use the experimental lattice constants a = 3.91 A and c = 24.984 A. Since the internal 
parameters of this structure have not been determined experimentally, we evaluated the 
equilibrium positions of the atoms using a damped Newton dynamics method. (The final 
force on each atom was less than m 1 mRy/a.u.) The equilibrium coordinates are listed in 
Table I. 

The band structure of the optimized geometry is shown in Figure 2 along several high- 
symmetry lines in the Brillouin zone. Roughly speaking, the two bands centered at ~ — 10 
eV are derived from Bi 6s orbitals while the manifold of 27 valence bands are derived mainly 
from O 2p orbitals. We found that the fundamental band gap is indirect, since the valence 
band maximun lies at X while the conduction band minimun is the Brillouin zone center 
T. As is typical in LDA calculations for semiconductors, the band gap is underestimated. 



The experimental band gap obtained from UV-absorbance measurements is 4.2±0.2 eV |20 
which is twice as large than our theoretical value of ~ 2 eV. It is interesting to note that, in 
disagreement with our results, an indirect (T-M) band gap was obtained by the tight-binding 
calculation ||12|| , which indicates important differences in the details of the dispersion. 

The total density of states between -35 and 5 eV is shown in Figure 3 illustrating the 
energy position of the semicore states. The calculated total DOS is in reasonable agreement 
with XPS data, which show a valence band width of m 7 eV and a peak centered ~ 10 eV 
below the valence band maximum (Bi 6s band) 
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Although the valence band presents mainly 2p character, there is a quite strong hy- 
bridization with Bi and Ta states, as is evident from the sphere-projected density of states 
(DOS) shown in Figure 4. Examination of the DOS reveals that there is substantial O 
2p character in the conduction bands, rising from zero at the conduction band minimun 
with increasing energy. Conversely, there are strong Ta and Bi contributions to the valence 
band. The Ta 5d contribution is zero at the valence band maximun but rises strongly with 
decreasing binding energy, reflecting the Ta 5d-0 2p covalency (as it is common in ABO3 
perovskites). The Bi contribution has both s-like (mainly in the upper part of the valence 
band) and p-like (mainly in the lower part) character. These contributions arise from a 
strong hybridization between 0(2,3) 2p valence band states with Bi 6s (fully occupied) and 
Bi 6p (conduction) bands. A detailed examination of the band characters reveals that the 
valence band maximun and the conduction band minimun are not mainly localized in the 
Bi-0(3) layer, as was obtained by the tight-binding calculation |T^]. In fact, in our case the 



valence band maximun at X is primarily of 0(1,4) 2p character. 

We finally show in Figure 5 the valence charge density map of SBT in two high symmetry 
planes, where the Ta-0 and Bi-0 "covalent bonds" can be seen. 

B. Ferroelectric instability 

The main features of SBT described so far resembles the electronic structure of PbTiOs. 
In this material, besides the large Ti 3d-0 2p hybridization, there is an important hy- 
bridization between Pb 6s and O 2p states; and this covalent bond plays a central role in 
the stabilization of the tetragonal ferroelectric structure of PbTiOa |^T] . 



The presence of quite strong Ta-0 and Bi-0 hybridizations in SBT opens fundamental 
questions about the origin of its ferroelectricity: What is the underlying static potential for 
the ferroelectric distortion? Which is the primary source for ferroelectricity? We investigate 
these questions by a lattice dynamics study using the frozen-phonon approach. 

To search for the presence of a possible lattice instability of the tetragonal structure, we 
determined the phonon frequencies and eigenvectors of the infrared-active modes, which 
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are polarized perpendicular to the c axis. To this end, we calculated atomic forces for several 
small displacements (~ 0.01 A) consistent with the symmetry of the mode. From the force 
as a function of displacement, the dynamical matrix was constructed and diagonalized. The 
calculated frequencies and eigenvectors are listed in Table II. 

While experimental studies were carried out on infrared-active (IR) phonons in or- 
thorhombic SBT, using reflectivity and transmission measurements ||22|| , there is no ex- 



perimental determination of IR phonon frequencies for the tetragonal structure to directly 
compare with. The remarkable fact of our calculation is, however, the presence of one unsta- 
ble phonon mode. Roughly speaking, this mode mainly involves movements of the Bi atoms 
with respect to rest of the lattice. Actually, its displacement pattern is (0.085, 0.114, -0.206, 
0.325, 0.549, 0.053, 0.224, 0.182), showing large displacements of the TaOg perovskite-like 
blocks relative to Bi atoms. This displacement vector is obtained from the eigenvector by 
dividing each component by the square root of the corresponding atomic mass and then 
normalizing to unity. So, the polarity of the mode can be described as a vibration of the 
TaOg perovskite-like block relative to Bi, with an additional contribution arising from the 
movement of Ta relative to its surrounding oxygens. Rae et al. H argued that the polarity of 
orthorhombic SBT can be quite well described as a movement of Ta, 0(4) and 0(5) relative 
to Bi and 0(3) of the Bi02 layer with and additional movement of Ta relative to 0(4) and 
0(5). Although the displacement pattern of the unstable mode is not directly related with 
the atomic positions of the orthorhombic phase, both pictures show some consistence. It is 
worth to mention that, from a detailed analysis of the calculated force constants, the lattice 
instability primarily arises from the attractive Bi-0(2) interaction. 

The resulting in-phase movement of the Ta atoms with respect to their oxygen octahedra 
(obtained from the eigenvector of the ferroelectric mode), could be the explanation of the 
low remanent polarization observed in SBT. However, this situation could be different for 
the case of SBN, considering the different structural behavior of KNbOs and KTaOs. While 
KNbOs has a series of ferroelectric phase transitions, KTaOs remains cubic down to low 



temperatures. Furthermore, it was emphasized by Singh [p3l , that the absence of ferroelec 



tricity in KTaOs is due to the extreme sensitivy of the soft-mode to the covalency and the 



slight chemical differences of Nb and Ta, particularly the higher d binding energy of Nb. 
This point will be clarified in the future by a corresponding investigation on SBN. 

Finally, the total energy is evaluated as a function of the relative displacement of Bi 
with respect to 0(2) corresponding to the unstable mode. The results are shown in Figure 6 
which depicts the energy per formula unit for displacement patterns along the [100] and [110] 
direction, respectively. Ferroelectric instabilities with energy gains of ~ 4 mRy/cell and 6 
mRy/cell for the [100] and [110] directions are observed. Although the ferroelectric mode 
mainly involves displacements of the Bi atoms with respect to rest of the lattice, a [110] 
displacement of the Bi sublattice alone does not produce a lattice instability (see Figure 6). 
This indicates that Bi has not a tendency to go off-center in the tetragonal phase, and the 
energy wells presented in the figure are indeed associated to the specific pattern of atomic 
displacements of the ferroelectric mode. 

As already mentioned, the bonding of these materials stems from the transition metal- 
oxygen hybridization. As can be seen from Figure 7, the whole spectrum of eigenvalues gets 
modified in an intricate way upon the complex distortion that gives rise to the instability. 
The two main effects are an increase in the band gap and a decrease of the separation 
between the Bi-s states and the valence band. When the distortion sets in, spectral weight 
from the lower part of the conduction band gets reduced and transfered to higher energies. 
This is more pronounced in the Bi contribution to the DOS than in that of Ta. In the valence 
region, the DOS gets reduced at the top and increased at the bottom. Here, the effect on 
the Ta contribution is stronger. So, both covalent bonds seem to play an important role. 
It is worth to mention however that the energetics shown in Figure 6 finally comes from a 
very delicate balance between several contributions: covalent, Coulomb and repulsive ionic 
interactions. 

The [110] displacement results in a deeper total energy minimun, which is consistent with 
the fact that the low temperature ferroelectric state of SBT has orthorhombic symmetry. 
The existence of a saddle point in the [100] direction on the total energy surface could 
indicate that the phase transition in SBT is not as simple as considered until now. In 
order to discuss this point it is useful to remark some aspects of the dynamical mechanism 

8 



leading to the sequence of phase transitions in ABO3 perovskites. Theses crystals have 
been considered for a long time as displacive-type ferroelectrics. The main evidence for this 
behavior has been the existence of a T-TO soft mode which has been observed in many 
perovskites P^ES]. However, there is experimental [EB-CT as well as theoretical [|5T|^ |U3[ 



evidence that the sucessive phase transitions in KNbC>3, and BaTiC>3, have pronounced 
order-disorder features, and that they are quite well described in the framework of the 
eight-site model. According to this order- disorder approach, the total energy surface has 
a maximum for the cubic perovskite structure, eight degenerate minima for the [111] soft 
mode amplitude displacements, and saddle points for the [100] and [110] displacements. In 
the cubic phase, the eight minima are occupied with equal probability, where this symmetry 
is broken as the temperature is lowered: Four sites are occupied in the tetragonal phase, two 
sites in the orthorhombic phase, and finally, only one site is occupied in the rhombohedral 
structure. In this way, the relaxation process which displays a critical slowing down when 
T c is approached in the different phases is interpreted to be the driving mechanism of the 
phase transitions. 

The lattice instability related to the ferroelectric phase transition of SBT was recently 
studied by Raman scattering in the temperature range between 293K and 958K ||34j| . The 



lowest frequency mode of 29 cm -1 at room temperature showed remarkable temperature 
variations towards T c . While its frequency decreases markedly, the damping factor increases 
rapidly below T c . This fact may suggest that the nature of the phase transition shows 
a crossing over from displacive to order-disorder type in the neighbourhood of T c . As in 
KNbC>3, it was also found, that the extrapolated frequency is still finite at T c , where this 
relatively large value may originate not only from the first order phase transition, but also 
due to the coupling between the soft mode and the strain. 

If this coupling between soft mode and strain is strong enough, the energetics shown in 
Figure 6 would give rise to the presence of an intermediate phase in SBT, leading to the 
following phase transition sequence: paraelectric (tetragonal) — > ferroelectric (net polariza- 
tion along the [100] direction — > ferroelectric (orthorhombic, net polarization along the [110] 
direction). As in the eight-site model, the four energy minima would be occupied with equal 
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probability in the high-temperature paraelectric phase, two sites would be occupied in the 
intermediate phase (with a net polarization along the [100] direction), and finally, one site 
would be occupied in the orthorhombic phase. 

Recently, an anomaly at T=520K was observed in specific heat measurements of SBT 
films ||35|| . For Bi-rich Sro.8Bi2.2Ta20g two anomalies were found at 620K and 410K, with 
corresponding changes in the X-ray diffraction pattern suggesting structural changes. The 
introduction of a small amount of excess Bi improved significantly the ferroelectric properties 
of SBT (the spontaneous polarization is two times larger and the Curie temperature shift 
to 670K from 608K ||36|| ). These results also suggest that the phase transition in SBT is 



not as simple as considered until now. More detailed examinations of crystal symmetry 
and dielectric measurements on single crystals are necessary for clarifying the nature of 
ferroelectricity in the SBT family compounds. 
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Figure Captions 



Figure 1: 

Tetragonal structure of SrBi 2 Ta 2 9 . Only atoms between |c and |c are shown. 
Figure 2: 

Band structure of SrBi 2 Ta 2 9 along several high-symmetry lines in the Brillouin zone. 
Figure 3: 

Total density of states for SrBi 2 Ta 2 9 , between -35 and 5 eV, showing the energy position 
of the semicore states. 

Figure 4: 

Total and site-projected electronic densities of states for the valence and conduction bands 
of SrBi 2 Ta 2 9 . 

Figure 5: 

Valence charge density plot in the (100) and (110) planes of SrBi 2 Ta 2 9 . The scale is from 
to 0.2 electrons per bohr 3 and the contour interval is 0.01 electrons/bohr 3 . 

Figure 6: 

Energy as a function of the distortion along the [110] (•) and [100] (□) directions corre- 
sponding to the unstable mode. The normal coordinate is represented by a motion of Bi 
relative to 0(2). The energies are with respect to that of the perfect tetragonal structure. 
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The energy as a function of the [110]-displacement of the Bi sublattice is also shown (A). 
Figure 7: 

Comparison of the partial density of states for Ta and Bi in the tetragonal structure, with 
(full line) and without (dotted line) displacements of the ions, according to the displacement 
pattern of the unstable- mode along the [110] direction. 
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TABLES 



TABLE I. Equilibrium atomic coordinates in lattice constant units (a=3.9lA and c=24.984 A). 



atom 


A 


Y 


L 


br 


o.o 


0.0 


0.0 


Ta 


0.5 


0.5 


± 08468 


Bi 


0.0 


0.0 


± 0.20124 


0(1) 


0.5 


0.5 


0.0 


0(2) 


0.5 


0.5 


± 0.16033 


0(3) 


0.0 (0.5) 


0.5 (0.0) 


0.25 


0(4) 


0.0 


0.5 


± 0.07602 


0(5) 


0.5 


0.0 


± 0.07602 
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TABLE II. Frequencies oj and eigenvectors of the modes, which are polarized perpendicular 
to the c axis. 



UJ 








Eig 


envector 








(cm -1 ) 


Sr 


Ta+ 


Bi+ 


01 


02 + 


03+ 


04+ 


05+ 


495 


0.036 


0.052 


0.001 


0.432 


0.033 


-0.001 


-0.617 


0.148 


247 


-0.051 


0.064 


0.053 


-0.754 


0.324 


-0.205 


-0.212 


0.125 


215 


-0.032 


0.053 


-0.117 


-0.360 


-0.122 


0.621 


-0.113 


0.073 


113 


-0.786 


0.395 


-0.026 


0.019 


-0.120 


-0.066 


-0.019 


-0.122 


100 


0.147 


0.037 


0.031 


0.010 


0.239 


0.108 


-0.129 


-0.634 


73 


0.502 


0.301 


-0.191 


-0.232 


-0.403 


-0.199 


-0.088 


-0.103 


i54 


0.130 


0.251 


-0.489 


0.213 


0.360 


0.035 


0.147 


0.120 
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